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Abstract 

_C A survey on the dynamical and thermodynamical properties of plasmas with 

strong Coulomb interactions in the quasi-classical density-temperature region is 
given. First the basic theoretical concepts describing nonideality are discussed. 
The chemical picture is introduced. It is shown that the nonideal plasma sub- 
On' system of the free charges has a rather large quasi-classical regime, where the 
quantum effects yield only corrections to the merely classical dynamics. The 
plasma of free charges may be described by effective potentials which incorporate 
quantum effects in an approximative way. The simplest effective potentials are 
Qh| only space-dependent, more advanced methods include momentum-dependent 
interactions. On the basis of these potentials analytical results are derived and 
simulation methods are developed. It is shown that effective potentials are ap- 
propriate for the description of thermodynamical as well as collective properties. 
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1 Introduction 

and in technology |], ||, |3|, || . In these plasmas the mean potential energy is of the same 
order of magnitude as the mean kinetic energy. Then we speak also about nonideal 
plasmas. We will study in this work the dynamics, thermodynamic properties and 
several collective effects of strongly coupled non-degenerate one-component plasmas 
(OCP) and symmetrical two-component plasmas (TCP). The investigation is restricted 
to the subsystem of the free charges, which is defined by means of the chemical picture 
0. In this model the atoms, ions, and molecules are treated as separate species. 
Therefore the constituents of the plasma are free electrons, free nuclei, ions, atoms, and 
molecules. All species are treated on equal footing (principle of particle democracy). 
*email:werner@summa.physik.hu-berlin.dc 



The advantage of the chemical picture is that it is in many cases more appropriate for 
the description of real plasmas . In the chemical picture the free charges constitute 
a well defined subsystem of the plasma. 

The main purpose of this work is the investigation of a quasi-classical dynamics 
of free charges based on effective potentials which may be used to describe thermo- 
dynamical and collective properties of the free charges in this subsystem. We restrict 
ourselves to the classical and near-classical region. Degeneracy is taken into account 
only in an approximative way. Therefore the new results refer in particular to a region 
where the plasma is still nondegenerated but nevertheless strongly coupled 0. In 
this work we will restrict ourselves to one-component plasmas and to two-component 
plasmas which are anti-symmetrical with respect to the charges (e_ = — e+) and sym- 
metrical with respect to the densities (n+ = rii = n_ = n e ). In particular we consider 
the model case of electron - positron plasmas m + = m e and H-plasmas with the mass- 
relation m + = 1840 m e . We include the - so far unrealistic - case of mass - symmetrical 
plasmas since in this case the thermodynamic functions and other analytic functions 
describing the plasma are of particular simplicity. This is due to cancellation effects 
caused by the symmetry of the masses and anti-symmetry of the charges 0. The 
effective potentials may be used to calculate correlation functions, thermodynamic 
properties and structure factors of the free charges in semi-classical non-degenerate 
quantum plasmas. The effective potentials are first obtained from the Slater sum 
method. Then momentum-dependent potentials are introduced and discussed. 

Characteristic parameters: The average distance of the electrons is the Wigner- 
Seitz radius d = [3/47rn e ] 1 / 3 , the Bohr radius is defined as a# = h/me 2 . Other 
characteristic lengths are the Landau length I = e 2 /kT, the De Broglie wave-length 
Aj = h/^TTrriikT] 1 ' 2 of particles of species i and thermal wave length of relative motion 
Ay = h/(2mijkT). Furthermore we define the following dimensionless parameters: 
(i) the coupling strength: T = l/d. 

(ii) the degeneration parameter: n^A 3 or 6 = 2mkT/[h 2 (37r 2 n e ) 2 / 3 ]. 
(iii) the interaction parameter £ y - = — (eiCj) / (kT \j) . 

Let us first discuss the OCP: In the limit when the average distance is large in 
comparison with the Bohr radius, i.e. if d > a% the electron gas behaves classical. The 
classical case was treated analytically by Abe and others || , Monte Carlo calculations 
in a wide range of T-values were carried out e.g. by Brush, Sahlin and Teller, DeWitt, 
Ichimaru and other workers (see |], |], [| |J). Quantum corrections to the classical 
case which are relevant at moderate values of r s were investigated by many authors 
fl], |H], I!] . All the analytical calculations mentioned so far cover only limiting cases 



as e.g. r s < 1 or T < 1. This is the reason why simulations are of large interest. 
There exists extensive Monte Carlo (MC) calculations for the classical region [12 



Classical molecular dynamics (MD) calculations were presented by Hansen et al. |13[ 



The particular interest in MD calculations is connected with the fact, that they give 



us also access to non-equilibrium properties [|TJ[ . 

For the TCP the available material is less exhaustive. We mention analytical calcu- 
lations in the region of small T |5|, |6|, [15], [L6| and Pade- approximations |2|, [T7[] connecting 



the analytically accessible regions. Further we mention several investigations devoted 



to the simulation of two-component plasmas [|18| , |19| , |20| , |2l] , |25| , P8[ . In particular our 
interest is devoted here to quasi-classical methods. Quasi-classical simulations of two- 
component plasmas were pioneered by Norman and Hansen [18|, [H|. These methods 
attracted a great deal of interest because of their relative simplicity. A special di- 
rection in developing simulation methods uses momentum- dependent potentials. The 
main idea in this approach is to model quantum effects by certain constraints in the 
phase space constructing an appropriate Hamiltonian |22|, ^3], [24]]. In particular the 
Pauli exclusion principle is simulated by a momentum- dependent two-body interac- 
tion. 

2 Effective Interaction Potentials 

2.1 Space - Dependent Effective Interactions 

As pointed out the idea of quasi-classical methods is to incorporate quantum-mechanical 
effects (in particular the Heisenberg and the Pauli principles) by appropriate poten- 
tials. Such a quasi-classical approach has of course several limits which are basically 
connected with the trajectory concept. We mention for example the principal difficulty 
to describe microspic quantum effects as tunnelling, and macroscopic quantum effects 
as superfluidity and superconductance. Our aim is only the calculation of standard 
macroscopic properties which have a well defined classical limit. Since bound states 
cannot be described classically our methods are restricted to the subsystem of free 
charges. However, this is not a very serious restriction since most of the plasma prop- 
erties are determined by the subsystem of the free charges. 

The easiest way to arrive at effective potentials describing quantum effects is the use 
of the so-called Slater sums which are defined by the N - particle wave functions, 

5(ri,...,rAr) = const ^exp(-/5E„) |* n (ri,...,rAr)| 2 . (1) 

The integrals over the distributions S(ri, . . . , r?v) yield the correct quantum sta- 
tistical partition function. The Slater sums for Coulombic systems were studied in 
detail by several authors |j], Q . With the knowledge of the Slater sums one gets exact 
space distributions in equilibrium by the choice 

U^ N \r 1 ,...,r N ) = -kT\nS(r 1 ,...,r N ) . (2) 

These potentials are often called quantum statistical effective potentials and they 



are used to calculate the partition function pi 0, M. Thus one obtains the correct 



thermodynamic functions. However, a severe disadvantage of this choice of the effective 
potentials eq.(Q) is, that the momentum distributions are correctly described only in 
the Boltzmann limit. Therefore this approach is principally restricted to the non- 
degenerate region. 

The Slater sum is an analogue of the classical Boltzmann factor and one defines 
therefore an effective potential by 



S%!(r) 



exp {-(3u ab (r)) = const. J^exp {-(3E a ) \ ^ a 



(3) 



Here \I/ Q and E a denote the wave functions and energy levels of the pair ab. 

A quantum mechanical calculation including the first orders in the perturbation 
theory was first given by Kelbg; a similar more simple effective potential was derived 
by Deutsch and was used in the simulations by Hansen and McDonald Jl9| . 
The effective potentials derived from perturbation theory do not include bound state 
effects. In order to treat the region, where bound states are of importance, a quite 
different approach is necessary |J. The first step is a transition to the chemical picture 
i.e. bound and free states have to be separated. The second step is the derivation of 
an effective potential for the interaction of the free charges in the plasma. In order to 
proceed we split the Slater sum into a bound part and a free part. This splitting is 
not unique. We use here the so - called Brillouin - Planck - Larkin (BPL) convention. 
This way of division free - bound corresponds to a smooth cut - off of the partition 
function which is based on omitting of the divergent elements in the sum. This leads 

to! 



S b ab (r) = const.J2' (exp (-(3E a ) -l + /3E a 



9, 



(4) 



where the sum extends only over all the discrete states. It is well known, that the 
BPL - convention leads to the simplest expressions for the thermodynamic functions. 
The Slater sum of the free charges is defined by [p], §, §]: 

S: b (r) = S$-S b ab = exp(-f3u: b (r)) . (5) 

In this way we arrive at an effective interaction potential of the free particles u* b (r), 
which is finite for h ^ and has a weak (integrable) singularity in the limit H — > 0. 

For the electron - electron and for the ion - ion interaction the classical limit gives 
simply the Boltzmann factors, i.e. u* ee = u* { = e 2 jr. For the electron - ion pairs one 
has to perform the classical limit in eq.(|5|). Explicitly one finds for the classical limit 
within the BPL - convention f|6l fl : 
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(6) 



where $(x) denotes the error function. The effective potential 

u: b (r) = -kT-log(S: b (r)) (7) 

may be used for simulations of the free charges in the purely classical region. 



2.2 Momentum - Dependent Effective Interactions 

As mentioned already above, a principal disadvantage of purely space-dependent po- 
tentials is the incorrect representation of the momentum - distributions of the plasma. 
In order to achieve a correct representation of the Fermi distribution for the momenta, 



momentum-dependent potentials have to be included [^J, |2^, |27|. In what follows we 



will assume a quasi-classical Hamiltonian of the following structure: 

H = Jl^ + Y.Vp(— >— )+Y.--f{— >— ) • (8) 

l=1 2m t<j \r ij0 PijoJ i<3 r \r ij0 p ij0 J 

Here r^- is the usual distance in the coordinate space and pij the distance in the 
momentum space. Further we define the characteristic parameters 

2 2,2 2 2,2 f- / /n\ 

r ijo = r io + r jo; Pijo = Pio + Pjo', Pio = h/no , (9) 

where r^ is a characteristic length (i.e., the radius of the wavepacket) of the par- 
ticle i. We have in the Hamiltonian two kinds of particle interaction: the so-called 
Pauli-potential Vp acting only between identical particles and a Coulomb interaction 
modified by a certain function F(x,y). In order to derive effective expressions of this 
type the Hamilton operator H is averaged with respect to test wave functions E 



H(q,p;h) = j ' dxil)* (x)Hil> {x) . (10) 

This definition of an effective Hamiltonian stems from the so-called wave-packet dy- 
namics. In the last time it has found several applications to plasmas |25], |26|, |29|, |3(J . 
If one chooses for the test wave functions symmetrized and anti-symmetrized combi- 
nations of minimum uncertainty wave packets for the particles of species k, 

ipko(x) = const exp f —^ h — J , (11) 

one does not get any Pauli potential for two electrons with antiparallel spins, for 
two electrons with parallel spins the following Pauli potential is obtained, 

V Pij (x,y) = 6 ij ^ T - exp(-A 2 ) A , . (12) 

mrfj v ' 1 — exp(— A z ) 

This is a two-body interaction depending on the phase-space distance A 2 = x 2 +y 2 . 
In our simulations we averaged over the two spin configurations and used a simplified 



Pauli-Potential which is purely Gaussian [53 



h 2 

V P (p,r) = V -exp(-A 2 ) , V ij0 = - . (13) 

zm r ^ 

The averaging over wave functions leads also to a modified Coulomb interactions, 
described by the function F (eq. |). We obtain 

F(x,y) = erf (x) . (14) 



This potential was obtained already by Klakow et al. ||25| , it represents the elec- 
trostatic energy between two charges which are Gauss-distributed. The free potential 
parameter r$ was fitted in such a way, that the properties of an electron gas without 
interactions as the binary correlation function and the Fermi momentum distribution 
are well reproduced. We assume that the momentum uncertainty is given by that of 
the free Fermi gas. In this way we get 

Pi0 = 3™4t ( 15 ) 

and Tjo = h/pio, where e^ n is the mean kinetic energy of the free Fermi gas of 
particles of species i. These relations yield at high temperatures p 2 = rriikT; r ii0 = 
2rf = h 2 /2rriikT = A^/2; Vq = kT/2. We mention that there exist other estimates for 
the free parameters of the Pauli-potential p3|. 



3 Quantumstatistical Theory and Simulations 
3.1 Thermodynamic Properties 



Quantum corrections to the classical electron gas were derived in earlier work HlO| [27 . 
Generalizing the methods developed earlier [^7| to the case of the TCP we assume 
that the interaction part of the free energy density of the plasma can be split into a 
classical and a quantum-mechanical part 

fint = fcl + fqu , (16) 

where f c i is the known free energy density of the classical plasma of free charges 
and f qu is the difference between the full and the classical free energy density. Explicit 
calculations for the classical free energy of the electron gas were given first by Abe 
for the low density case and extended by Cohen and Murphy ||. For a symmetrical 
classical TCP it was shown that (with the BPL convention) the first correction beyond 
the Debye term in the density expansion of the free energy vanishes. 0. 
Consider now the quantum-mechanical corrections to the classical free energy density. 
For the plasma of free charges this expression is convergent. In the low density limit 
we get 

f gu = -kTj2 n i n j$ B ij( T ) ■ ( 17 ) 



Here SB^ (T) is the difference between the quantum and the classical second virial 
coefficients. This function can be calculated exactly by using methods of quantum 
scattering theory |I5| . First the two particle trace is transformed to a contour integral 



over the resolvent of the Coulomb scattering problem. For a mass-symmetrical TCP 
all terms except to one cancel and we get (neglecting exponentially small degeneracy 
effects) the rather simple result 

f qu = -n^ ni n e \Ul + 0(n^) (18) 

If the masses are not symmetrical, infinite series in the £ - parameter appear, the so - 
called quantum virial functions Q and E |IL |J . The theory may be extended to higher 
concentration by using the method of Pade - approximations flTj. 



We carried out extensive simulations for the electron gas |27|] an d also a few simu- 
lations for the symmetrical TCP by using the effective Hamiltonian described above. 
Since our modified Coulomb interaction differs from the bare Coulomb interaction only 
for short distances the Ewald sum technique for handling the long range part could 
be used. To study equilibrium properties, as e.g. the average energy we carried out 
several Monte Carlo runs of 2 • 10 6 steps for ensembles of 64-512 particles. The results 
were extrapolated to infinite particle numbers. The regions of degenerate and nonde- 
generate plasmas were investigated. The results of our simulations were discussed in 



detail in earlier work [27, 2q. As the main result we may quote, that at all densities 



the deviations from the classical calculations are very small [pqj . Further we may 
state that the overall agreement between the analytical formulae and the simulations 
is rather good. The deviations are all within the error bars of the simulations. Only 
at low temperature a systematic deviation is observed, here the semi - classical model 
does not work. At conditions of weak or moderate degeneracy our model yields quite 
reasonable results. 

3.2 Collective Excitations 

In order to study collective nonequilibrium effects we simulated a system of 250 elec- 
trons with periodic boundary conditions by Molecular Dynamics (MD) calculations. 
Usually MD simulations consist of two parts: First the desired temperature is ad- 
justed using some thermostat, while in a second phase the energy is kept constant 
and measurements are performed. Since most known thermostats do not work in the 
case of momentum dependent interactions we replaced the equilibration phase by a 
MC simulation using the Metropolis Algorithm. This Algorithm is independent from 
any particular form of the interaction. To perform the MD simulations we used a 4th 
order Runge-Kutta integrator with stepsize control. The runs were of length of order 
1000 times the inverse plasma frequency. 

First the individual motion of the electrons was studied by calculating the velocity 



autocorrelation function < v (t + r)v(t) > t |27|, [30]. It was shown that for moderate 



coupling (r = 1) the velocity autocorrelation falls monotonically to zero , whereas 
for the case of strong coupling (T = 100) the velocity acf shows oscillations with a 
frequency close to the plasma frequency. 

To describe the collective behavior of the system we have investigated the dynamic 
structure factor of the electron system, 

S(k, u) = — j /_ oo e wt < p(k, t) p(-k, 0) > dt , (19) 

where p(k,t) = J2i ex P{~ikfi) is the Fourier component of the microscopic electron 
density. We obtained the dynamic structure factor from the MD simulations by ap- 
proximating the Heisenberg operator rl(t) by the position of the i-th particle in the 
simulations. However, the thus obtained quantity (we denote it by R(k, u)) is symmet- 
ric with respect to the frequency. It corresponds therefore to a classical fluctuation- 
dissipation theorem, 

R$,w) = {mr<j>(k)/3u)- 1 Ime- 1 (ic,u) . (20) 

where e(k,u) is the dielectric function of the electron OCP. It can be seen from Eq. 
(p0|) that R(k, u) cannot be regarded as a structure factor, but as a normalized loss 
function. In what follows we will discuss the normalized loss function. Notice, that in 
the classical case the loss function and the dynamic structure factor coincide. 

We have plotted the loss function R(q, uS) (q = ka) for various q values, for the 
case of moderate coupling (r = 1), strong coupling (r = 10 ) and for the case of very 
strong coupling (r = 100 ) and for different parameters of degeneracy 9 = 1 (moderate 
degenerate) and 9 = 50 (classical plasma) (Figs. 0-0). For the regime of moderate 
coupling we have compared the results of the simulations with the corresponding data 
from the Random Phase approximation (RPA), |30[. We see from Figs. [I], |2] that 
in the case of a moderate coupled plasma (r = 1) the shape of the loss function 
calculated from the simulations is damped stronger and slightly shifted to the left in 
comparison with the RPA loss function. For moderate coupling and in both cases 
(weakly degenerate and classical plasma) the plasmon peak can be observed only for 
the smallest q value (q=0.619). From Figs [l] and ^| it can be also seen that the change 
of the degeneracy parameter 9 in the range from 50 to 1 has only a small influence on 
the results. 

However, at higher degrees of degeneracy (9 = 0.15) the plasmon peak obtained 
from the MD datas is shifted towards higher frequencies. But the peak position in this 
case differs quite significantly from those predicted by the RPA which is shifted much 
more due to the increase of the average velocity (the Fermi velocity) which leads to a 
strong positive dispersion (Fig. |3|). The simulations underestimate this shift. 

In the strong coupling regime (r = 10) and for a moderate degenerate plasma 
(9 = 1) we observed at the smallest q value a sharp plasmon peak centered near up 
(Fig. |]). With the increasing wavevector the plasmon peak widens, a plasmon peak 



can be observed up to q = 1.856. At still higher q values the peak vanishes. Notice 
also that within the resolution limited by the simulations almost no dispersion could 
be observed at this thermodynamic conditions. We interpretate this as the limiting 
region from positive (at smaller T) to negative dispersion (at higher T). 

In the case of very strong coupling T = 100 an extremely sharp plasmon peak 
centered near ojp can be observed at the smallest q value (Figs. |^, ||). Here the 
plasmon peak can be observed up to q = 3.094. For both the case of a classical plasma 
(6 = 50) and for the case of moderate degenerate plasma (6 = 1) a negative dispersion 
is seen (with increasing wavenumber the peak position is shifted more and more to the 
left). This behavior contradicts to the RPA where no plasmon peak is predicted in 
this regime due to the strong Landau damping . However, the RPA cannot be applied 
to the strong coupling regime. On the contrary, the results of our simulations for the 
case of a weakly degenerate plasma (8 = 50) are in a good agreement with the results 
of corresponding MD simulations of Hansen et al. for the classical one component 
plasma [|l^]. From Figs. || and || one also sees that the negative dispersion is more 
pronounced in the case of a classical plasma. 

Further we have compared the results of our simulations with the expression of 
the dynamic structure factor obtained by the application of the classical theory of 
moments which is appropriate also for the case of strong coupling |3l| . A more detailed 



discussion can be found elsewhere | 3D| . Here we show the results for two different q 



vectors at T — 100 (Figs.[7|, || p0| ). The agreement of the loss functions calculated 
by the application of the theory of moments (sum rules approach) with that from the 
MD calculations is rather good. The theoretical curves reproduce the varying shape of 
the dynamic structure factor and describe the plasmon peak position in a satisfactory 
manner. However, the agreement in the height of the peaks is less satisfactory . 

Thus we can conclude that due to the quasi-classical character our quantum molec- 
ular dynamic simulations describe the collective excitations of the electron gas only 
approximately. Our model yields quite reasonable results at weak and moderate de- 
generacy, whereas for the case of high degeneracy it seems to break down. 

4 Discussion 

We developed here a simple quasi-classical model of quantum plasmas based on a dy- 
namics with an effective momentum-dependent Hamiltonian. The quantum-mechanical 
effects corresponding to the Pauli and the Heisenberg principle were modeled by con- 
straints in the Hamiltonian. By using the concept of minimum uncertainty wave pack- 
ets, momentum-dependent effective potentials were derived. The effective potentials 
were used to simulate one-component plasmas and mass -symmetrical two - compo- 
nent plasmas by means of MC and MD methods. The result of the simulations is in 
good agreement with analytical calculations of the thermodynamic properties in the 
region of small degeneracy and moderate coupling. 

9 



MD studies provide also informations on the dynamical properties of the plasma, 
in particular on collective excitations. As a basis to study these effects we calculated 
the structure factor of the plasma. The agreement between the simulations and the 
analytical theory based on RPA, sum rules and the theory of moments is reasonable. 
In most cases the shape and the location of the plasmon peak is reproduced in a rea- 
sonable way, in some cases larger deviations are observed. We cannot expect however, 
that the present model, which was obtained by fitting the potential to equilibrium 
properties will describe all non-equilibrium properties in a quantitative way. Further 
improvements of the model might be unavoidable. We may hope however that at least 
near to equilibrium some realistic features are still reflected by the model. 

Acknowledgment: The authors thank V.M. Adamyan, Yu.M. Klimontovich, B. 
Militzer, V. Podlipchuk and I.M. Tkachenko for helpful discussions. 
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Figure 1: Comparison of the MD and RPA loss function R(q,u) versus frequency 



uj/uj p for different wavevectors q at V = 1 and 9 = 1, [pO . 
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Figure 2: same as Fig. |1|; for T = 1 and 9 = 50, 
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Figure 3: The MD loss function R(q, uj) versus frequency uj /uj p for wavevector q = 0.619 
at fixed T = 1 and different 8. 
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Figure 4: The MD loss function R(q, u) versus frequency u/u p for different wavevec- 
tors q at T = 10 and 8 = 1. 
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Figure 5: The MD loss function R(q,u) versus frequency uj/uj p for different wavevec- 
tors q at T = 100 and 6 = 50 
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Figure 6: same as Fig. [5], at T = 100 and 6 = 1. 
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Figure 7: comparison of the MD loss function R(q, to) versus frequency uj/uj p with 
the corresponding loss function from the sum rules approach at V = 100 and 9 = 50 
for wavevector q = 1.856, . 



4.0 



G O 



"quantum" MD results 
sum rules approach 



3.0 



o 
cr 



2.0 - 



cc 



9 & 



1.0 



0.0 




0.0 



0.5 1.0 1.5 

w ( in units of the plasmafrequency ) 



2.0 



Figure 8: same as Fig.|7|; at T = 100 and 9 = 50 for wavevector q = 3.094, . 



15 



